SUPPLEMENTAL MATERIAL 



TABLE I: parameters in G-protein cycle 



parameter 


value 


References and Notes 


kl 


0.12nAf mm"' 


2 X lO^Af-'s-' (12) 


k2 


0.6min~^ 


0.01 S-' (12) 


k3 


0.24mm-' 


4 X 10-^s-' (12) 


k4 


320nM ■ min 


nt to the time-course oi G protcm activation 
(12). (See Fig.2A ) 


k5 


7 nM 


k6 


30nM • min~^ 


4:{molecules/cell)s^ ^ (12) 


k7 


0.024min-' 


4 X lO-^'s"' (12) 


k8 


0.0048mm-' 


IQ-" {molecules /cellos'' (12) 


k9 


115nAf ■ min^^ 


fit to the time-course of G protein activation 
(12)(See Fig.2A), 

and experimental data that G protein increased from 
8000 molecules/cell to 12000 molecules/cell (8) 


klO 


60 nM 


kll 


20nM • mm-' 


kl2 


0.08mm-' 


kl3 


0.24mm-' 


4 X 10-^s-' (12) 


kl4 


0.025nM-'mm-' 


0.11s-' (12) 


kl5 


480nM-'min-' 


l(mo/ecit;e/ce//)-'s-'(12) 



TABLE II: parameters in recruitment of scaffolds 



parameter 


value 


References and Notes 


kl6 


0.05nA/-'mm-' 


ti/2 = 8.22s for SteSs' recovery 


kl7 


5mm-' 


at shmoo tips (49) 


kl8 


0.00007nA/-' min-' 


fit to (24) 


kl9 


0.001mm-' 


(See Fig.2C. ) 


k23 


30mm-' 


ti/2 ~ 2s for the recovery of ste5 in the nucleus. 
(49) 



TABLE III: dephosphorylation 



protein 


al 

(nAf-'mm-'') 


dl 

(min,-"') 


pl 

(min "') 


a2 

{nM~'^min^^) 


d2 

(min-"') 


p2 

(min "') 


Stell(MAPKKK) 


0.5 


24 


10 


1 


24 


10 


Ste7(MAPKK) 


0.5 


24 


10 


0.5 


24 


10 


Fus3(MAPK) 


0.2 


10 


30 


0.4 


20 


30 



2 



TABLE IV: phosphorylation in cytosol 



protein 


a3 


d3 


p3 


a4 


d4 


p4 




{nM~^min'^) 




(mm ■"") 






[min ^) 


Ste7(MAPKK) 


1 


36 


10 


1 


36 


10 


Fus3(MAPK) 


0.1 


36 


10 


0.1 


36 


10 



Notes: Equations of reaction in this part are all in the following form: S + E ^ SE ^ S' + E. 
As for the underset n: 

1: dephosphorylation of the once phosphorylated kinase; 2: dephosphorylation of the twice phosphorylated kinase; 3: once 
phosphorylation; 4: twice phosphorylation. Parameters are estimated in analogy to Ref.ll. 



TABLE V: Other reactions in the cytosol 



parameter 


value 


References 


k24 


1.2 nM'^min'^ 


ReL23 


k25 


24min"^ 


ReL23 


k26 


0.0015 nM^^min^' 


Ref.40 


k27 


0.1 min~^ 


ReL41 




FIG. 1: 27 kinds of scaffold-kinase complexes in solution and 27 kinds of scaffold-kinase complexes at the membrane 
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Bl + Stell ^B2 

B1 + Ste7pp ^ B5 

B2 + Ste7 ^ B8 

B2 + FusSpp ^ B15 

B3 + Fus3 ^ B13 

B4 + Stellpp ^ B9 

B5 + FusSpp ^ BIO 

B5 + Fus3pp <- B18 

Be + Ste7 ^ B16 

B7 + Stell ^ B15 

B8 + Fus3 ^ B20 

B9 + Fus3 ^ B21 

BIO + FusSpp <- B27 

B12 + Ste7 ^ B20 

B13 + Ste7pp <- S24 

B15 + 5te7 ^ B23 

B16 + 5tell ^ B20 

B17 + Stell ^B22 

B19 + Stellpp ^ B25 

B22 -» B26 

CI + Stell ^ C2 

CI + Ste7pp <- C5 

C2 -> C3 
C2 + Fus3 ^ C12 
C3 + Ste7pp <- CIO 
C4 + Siell ^ C8 
C4 + Fus3pp ^ C19 
C5 + Fus3 C17 
C6 + Stellpp ^ C13 
C7 + Stellpp ^ C14 
C7 + Ste7 ^ C19 
C8 + F«s3pp «- C23 
C9 + Fus3pp «- C25 
CIO + Fus3pp ^ C27 
C12 1- C13 
C13 + Ste7 ^ C21 
C14 + Ste7 ^ C25 
C15 + Ste7pp ^ C26 
C17 -» C18 
C18 + Stellpp — C27 
C19 + Stell ^C23 

C22 C24 
C24 -» C27 
Bl ;=±C1 
B4 ;=±C4 
B7 ^C7 
BIO ^ CIO 
B13 ^ C13 
B16 ^ C16 
B19 ^C19 
B22 ^ C22 
B25 ^ C25 



Bl + Stellpp ^ B3 

Bl + Fus3 ^ B6 

B2 + Ste7pp ^ Bll 

B3 + Ste7 ^ B9 

B3 + Fus3pp 1- B14 

B4 + Fus3 ^ B16 

B5 + Stell ^Cll 

B6 + Stell ^ B12 

B6 + Ste7pp ^ B17 

B7 + STe7pp ^ B18 

B8 + Fus3pp ^ B23 

B9 + Fus3pp 1- B25 

B11 + Fms3 ;=iB22 

B12 + Ste7pp ^ B22 

B14 + Ste7 ^ B25 

B15 + Ste7pp ^ B26 

B17 -» B18 

B18 + Stellpp ^B27 

B19 + Stell ^ B23 

B24 B27 

CI + Stellpp ^ C3 

CI + Fus3 ^ C6 

C2 + Ste7 ^ C8 

C2 + Fus3pp ^ C15 

C3 + Fus3 ^ C13 

C4 + Stellpp ^ C9 

C5 + Fus3pp ^ CIO 

C5 + i^us3pp 1- C18 

C6 + Ste7 ^ C16 

C7 + Stell ^ C15 

C8 C9 

C9 ^ CIO 

CIO ^ Cll 

Cll + Fus3 ^ C22 

C12 + Ste7 C20 

C13 + Ste7pp 1- C24 

C14 + Ste7pp ^ C27 

C16 + Stellpp ^ C21 

C17 + Stellpp ^C24 

CIS + stell ^ C26 

C20 ^ C21 

C22 -» C26 
C25 ^ C27 
B2 ^C2 

B5 ;=±C5 

B8 ^C8 

Bll ^Cll 

B14 R=; C14 

B17 ^ C17 

B20 C20 

B23 ^ C23 

B26 ^ C26 



Bl + Ste7 ^ B4 

Bl + F«s3pp ^ B7 

B2 + Fus3 ^ B12 

B3 + Ste7pp ^ BIO 

B4 + Stell B8 

B4 + FusSpp <- B19 

B5 + Fus3 ^ B17 

B6 + Stellpp ^ B13 

B7 + Stellpp ^ B14 

B7 + Ste7 ^ B19 

B9 BIO 

BIO + Fus3 ^ B24 

Bll + Fus3pp <- B26 

B13 + Ste7 ^ B21 

B14 + Ste7pp <- B27 

B16 + Stellpp ^ B21 

B17 + Stellpp ^ B24 

B18 + Stell ^ B26 

B21 -» B24 

B25 B27 

CI + Ste7 ^ C4 

CI + Fus3pp <- C7 

C2 + Ste7pp «- Cll 

C3 + Ste7 ^ C9 

C3 + Fus3pp <- C14 

C4 + Fus3 ^ C16 

C5 + Stell ^Cll 

C6 + Stell ^ C12 

C6 + Ste7pp <- C17 

C7 + STe7pp ^ CIS 

C8 + Fus3 ^ C20 

C9 + Fus3 ^ C21 

CIO + Fus3 ^ C24 

Cll + Fus3pp ^ C26 

C12 + Ste7pp ^ C22 

C14 ^ C15 

C15 + Ste7 ^ C23 

C16 + Stell ^ C20 

C17 + Stell ^C22 

C19 + Stellpp ^ C25 

C21 C24 

C23 ^ C25 
C26 C27 
B3^C3 

B6 ^C6 

B9 ^C9 

B12 ^C12 

B15 ^C15 

B18 ^ C18 

B21 ^C21 

B24 C24 

B27 C27 
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TABLE VI: parameters in MAPK cascade on scaffold 



Protein 


on 

(nM^^ 
min~^) 


off 


**on 
1 {nM-^ 
min~^) 


**off 


phospliorylation 
1 (mm~^) 


References and notes 


Stell 

(MAP 

KKK) 


1.2 


10 


1 


10 


20 


Stellpp can rebind to the scafold (27) 


Ste7 

(MAP 

KK) 


1.2 


0.1 




0.5; 
20 


46 


Ste7pp in the scaffold undergoes hyper- 
phosphorylation by FusSpp, which accel- 
erates Ste7's dissociation from the scaffold 


Fus3 
(MAPK) 


1.2 


10 




200 


250 


ti/2 = 0.3s for the recovery of Puspp in 
shmoo tips(10,57) 



TABLE VII: parameters in downstream reactions 



protein 


value 


References and Notes 


k28 


SnMmin"^ 




k29 


AOnM 




kSO 


O.Olmm"^ 




kSl 


0.5nMmin~^ 




k32 


lOnM 




k33 


2nMmin~^ 




k34 


7nM 




k35 


0.002mm" ^ 




k36 






k37 


lmin~^ 




k38 


10mm~^ 




k39 


lOmin"^ 




k40 


O.lnM^'-min'' 




k41 


O.Olmm"' 




k42 


O.OlnM-'mm"' 




k43 


O.lmm"'^ 




k44 


S2nMmin~^ 


fit the time course of Sst2 determined in experiment (8) 


k45 


InM 


the same as k44 


kl() 


OXmitun ' 


tlic sa.iii as k 1 1 


k47 




ii/2 ~ 4s for the recovery of FusSin nucleus (57) 
[Fus3,„]/[Fus3out] = 1.4 at i = (45) 


k48 




k49 




[FusSin] /[FusSout] = 2 after 60 min treatment of 
a-factor (45) 


k50 


15min~^ 



TABLE VIII: initial amount of proteins 





VdiU-t-i iiiv± 1 


n rt1"pci n Ti H tpt ptpti ppq 


Stc2 


1000 


8000 molc/ccll (8) 


Sst2 




9000 mnlp/rpll ("S^ 
^UUU liiUit; / Ceil I Oy 




1 nnn 

-LUUU 


snnn mnip/z-pii i'^?"^ 




1 nnn 

-LUUU 




CJfpl 1 


9nn 

zuu 


Ipqq tViQ-n 9000 Tnr\lp /r'pll n 


r LlbOO / 


9nn 

zuu 


±(JLcLi djlliULlllL Ui DLL./ bliUUiCl UC iChh Liii^lll ZUUU liiUiC / cilitl dUULlt y<J/0 Ul 


Stc7 


10 


J. tif50i7T, 




c»ri/~in^~ ^000 Tvi/~ilo /(^oll til n ri r~iOT/^T"P T\ ri o'l^Tn r\Y\ a 
dUULlL tJlJ'J'J iiiUlL. / l^L.ii I J- J ; cLiiU. UUiUi c JJlici UlliUllt; 

induction, [Fus'^in\/[Fus'^ou^ — 1-4 (45) 




300 


MAPKKK-P 


100 




MAPKK-P 


100 




MAPK-Pin 


100 




MAPK-Pout 


100 




Ste5in 


125 


(1) 


SteSout 


0.1 




Cdc28 


100 
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TABLE IX: mutants 



No. 


Mutant 


BcliaA-ior 


liiiploiiioutatiou 


Reference 


1 


Stei-'""^ (ciidocyto- 
impaired) 


G protein activation 

IpVpls MTl 

V Cio LI IJ 


k3=0.08 


Ref 12 


2 






kSI— kS2— k44— 


Ref 12 




heximide(synthesis 
inhibitor) 


closed down 


k4=k9=0 




3 


SST2A 


O LIJJCl oCllOlljl V i-ijy 

upon pheromone 
induction 


k44— f) [.'i''!/2l* n— n 
nM 


Ref 8 


4 


2 X SST2 


response upon 
pheromone reduces 


k44=64, 

[Sst2]t=o=500 nM 


Ref8 


5 


2 X G/37 


super-sensitivity to 
pheromone induction 


[G/37]t=o = lOOOnM 


Ref8 


g 


TiPO'fiti^T'p tppH nfiplrQ 




i,of:— — 

J\.ZiU JVZi 1 iVtJ i 






are cut down 


doesn't attenuate 

Wltll Lllllt: 


k44=0 




7 


1 cci^/^i Q "1" in rr 

n Vi o t: li n r V 1 i" p H 

SteTpp is slowed 
down 


cW^tiVcltHJli Ui X" U,oOJJJ|J 

witVi i"iTnp 






8 


dissociating speed 
for normal Ste7pp is 


activation of FusSpp 
decreases 


* * offxK = 20 






enVianppH 








9 


shuttle of the scaffold 


super-sensitivity to 


k22=k23=0 






Ste5 is cut down 


pheromone induction 






10 


Stell is continuously 
activated 


SteTpp is activated 
while PusSpp is re- 

TAfpQcpn HowTi 


[Stenpp]t=o 
200nM, [S'iell]t=o = 
\M APKKK — 

P]totoal — 

0, [a — factor] = 




11 


Ste7 is continuously 
activated 


FusSpp isn't 
activated 


[Ste7pp]t=o = 
210nM, [Ste7]t=o = 
0, [MAPKK 

P]totoal — 

0, [q — factor] = 


Ref 46 
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ODE functions: 

^t^^ = -kl[a - factor][Ste2\ + fe2[Ste2„et,„e] - fe7[Ste2] + + fe6 

d[Ste2a,cUve] ^ ^^j^ _ /actor][Sie2] - k2[Ste2aM^a] - fc3[S<e2„rfi„s] 

d[S^t2^ = fc44[5.el2,..,„.]'- _ ,4,[^^,2...„,] 
dt fc452 + [Stel2„rfi„e]2 

d[G„t] 



^ = -fe8[5te2„.«„.][G] + fel5[G„d][G^7l + nt'f rt'Tr'^ 12 " ^^^[G] + fell 
dt fclC + |Stel2actii.e] 



fe8[Sfe2„<.H„r][G] - fcl3[G„t] - fel4[Got] [5st2act>„c] 
dt 

^I^^ = fel3[G<:,t] + fel4[GQtl[Sst2„et,^el " fcl 5 [G,, d] [G,3 7] 
dt 

^t^^ = fe8[Ste2<.eti„e][G] - fel5[Gad][G^7] - fe40[G^7][Forlppout] + fe41[Farlpp„„tGs7]20 - fel8[G;37][Ste20] + fel9[G37Ste20] 
dt 

= -fel8[Gfj7][Ste20] + fel9[G;37Ste20] 

dt 

d[G07Ste2Ol 

^ '^^^ i = fel8[G,37][Ste20] - fel9[G07Ste2O] - fcl6[G07Ste2O]Bl + felTCl - fel6[G,37Ste20]S2 + fel7C2 - fel6[G07Ste2O]S3 + fel7C3 - 

fel6[G^7Ste20]B4 + fel7C4 - fel6[G^7Ste20]B5 + fcl7C5 - fel6[Gs7Ste20]B6 + fcl7G6 - fel6[G:37Ste20]B7 + fel7C7 - fel6[G^7Ste20]S8 + 
A;17C8 - fel6[G,-j7Sie20]B9 + fcl7C9 - fel6[G,37Ste20]B10 + A;17C10 - fel6[G,37Ste20]Bll + fcl7Cll - fel6[G,37Ste20]B12 + fcl7G12 - 
fel6[G,<,7Ste20]B13 + fcl7G13 - fel6[G07Ste2O]B14 + fel7C14 - fel6[Gfj7Ste20]S15 + fel7C15 - fel6[G37Ste20]S16 + fel7C16 - 
fel6[G,-n''Ste20]S17 + fel7C17 - fel6[G.<7Sie20]B18 + fel7C18 - fel6[G,37Ste20]S19 + fel7C19 - fel6[G,37S<e20]S20 + fel7C20 - 
fel6[G^j7Ste20]B21 + fcl7G21 - fcl6[G^i7Ste20]B22 + fel7G22 - fel6[G07Ste2O]B23 + fel7C23 - fel6[G07Ste2O]B24 + fel7C24 - 
fcl6[G^7Ste20]B25 + fel7C25 - fcl6[G37Ste20]B26 + fel7C26 - fel6[G07Ste2O]B27 + fel7C27 

''■^^^^^^^ = pljfifK[5tellpMAPifii'i<r-P]+o//ft:Kft:(C2+C8+Cll + C12+C15+C20+C22+C23+C26+S2+B8+Bll + B12 + B15 + B20 + B22+ 
dt 

B23+B26)-onjfifK[Stell](G4+C6+C7+Gl+C16+C17+C18+C19+C5+B4+S6+B7+Sl+B16+B17+B18+B19+B5)-fe26[Stell][_FMs3ppo„t] 

d[Stellp] ^ -alKKK[Ste\lp]{[M APKKK - P]o - [StellpMAPKKK - P] - [StellppMAPKKK - P]) + dlKKK[StellpMAPKKK - P] + 
dt 

p2KKK[StellppMAPKKK - P] 
d[StellpMAPKKK - P] _ 
dt ~ 

alKKK[Stellp]([MAPKKK - P]o - [StellpM AP K K K - P] - [StellppM AP K K K - P]) - (dlxKK + plKKK)[StellpMAPKKK - P] 

dfStellpp] ^ -a2KKK[Stellpp\{[MAPKKK - P]o - [StellpMAPKKK - P] - [StellppMAPKKK - P]) + d2KKK[StellppMAPKKK - 
dt 

P] - aS,KK[Stellpp][Ste7] + (d3KK +p3jcjc)[SfellppSte7] - aiKK[Stellpp\[Ste7p] + (d4jrK + p^k K)[StellppSte7p\ + ** offKKKiCi + CIO + 
C9 + C13 + C14 + C21 + C24 + C25 + C27 + B3 + BIO + B9 + S13 + B14 + B21 + S24 + S25 + B27) - * * onKKK[Stenpp]{Cl + C4 + C5 + 
C6 + C7 + G16 + C17 + G18 + C19 + Bl + B4 + B5 + B6 + B7 + B16 + B17 + B18 + B19) 
d[StenppMAPKKK - P] _ 

dt ~ 

a2KKK{Stellpp]{[MAPKKK - P]o - [StellpMAPKKK - P] - [StellppMAPKKK - P]) - (d2KKK + p2KKK)[StellppMAPKKK - P] 
= -a3KK[Ste7][Stellpp] + d3KK [Stellpp5te7] + plKK[Ste7pMAPKK - P] + offKK{C4, + C8 + C9 + C16 + C19 + C20 + C21 + 

dt 

C23 + C25 + B4 + B8 + B9 + B16 + B19 + B20 + B21 + B23 + B25) - onKK[Ste7](Cl + C2 + C3 + C6 + C7 + C12 + C13 + C14 + C15 + SI + 
B2 + B3 + B6 + B7 + B12 + B13 + B14 + B15) - fc24[Ste7] [F«s3o„t] + fe25[F«s3o„t5te7] 

d[Ste7gtellpp] ^ a3KK[Stellpp][Ste7] - (dSKK + p3KK)[StellppSte7] 
dt 

'^^^^^'^P^ = -uIkkGMAPKK - P]o - [Ste7pMAPKK - P] - [Ste7ppMAPKK - P])[Ste7p] + dlKA:[Ste7pMAP/!:/s: - P] + 

dt 

P?,kk[SI< llppSlil] - n4KK[Ste7p][Stellpp] + d4KK [SfellppSfe7p] + p2/f /f [Ste7ppMAPA:A: - P] 

d[SLc,pMAt A A - f] ^ ^^^^ [Ste7p]{[MAPKK - P]o - [Ste7pMAPKK - P] - [Ste7ppMAPKK - P]) - (dlKK +plKK)[Ste7pMAPKK - P] 
dt 

d[Ste7pStellpp] ^ a4^^[stellpp][Ste7p] - (d4jcjf + p4jf )[StellppSte7p] 
dt 

d[Ste7pp\ ^ -a2KK[Ste7pp\{[MAPKK - P]o - [Ste7pMAPKK - P] - [Ste7ppMAPKK - P]) + d2KK[Ste7ppMAPKK - P] + 
dt 

piK K[StellppSte7p\ - a3K [Ste7pp] [Ftts3o„t] + (d3.K + p3jf )[Ste7ppFMs3out] - a4K [Ste7pp] [FMs3pout] + (d4jf + p4a- ) [Ste7ppFMs3po,jt ] + * * 
o//jfjf (C5 + C10 + Cll + C17 + G22 + C24 + B5 + B10 + Bll + B17 + B22 + B24) + **o//^jf (C18 + G26 + C27 + B18 + B26 + B27)-fc27[Ste7pp] 

d[Ste7ppMAP_ff if - P] [Ste7pp] ( [MAPKK - P] o - [Ste7pMAPKK - P] - [Ste7ppMAPKK -P])-(d2KK +p2kk) [Ste7ppMAPKK - P] 

dt 

d[FMs3„^t] ^ _a3ft:[Ste7pp][F«s3o„t] + <i3ft:[Ste7ppFus3o„t] + pljf [Fus3po„«MAP/<- - P^ut] + o//ft:(C6 + C12 + C13 + C16 + C17 + C20 + 

dt 

C21 + C22 + C24 + S6 + B12 + B13 + B16 + B17 + B20 + B21 + B22 + B24) - ouk [FusSo^t] (CI + C2 + C3 + C4 + C5 + C8 + C9 + GIG + Gil + 

fe32rStel2 +■ 1^ 

Bl + B2 + B3 + B4 + B5 + B8 + B9 + BIO + Bll) - fe24[Ste7] [Fus3„^t] + k25[Fus3outSte7] + fe47[Fus3i„] - fe48[F«s3o„t] + — ^""^ 

fe5^ + [Stel2acUve] 

d[Fus3oMtStellpp] , „ , , 

— — = a3jf [S'te7pp][Ftts3o„t] - (d3if + p3/f ) [Ste7ppFtts3o„t] 

d[J »--.);;„„f] ^ ^alKlFus3po^t][MAPK - P„„t] + dlx [F«s3po„tMAPii: - P„„t] + p3jf [Ste7ppFus3o„t] - a4A[5te7pp][F«s3po„t] + 
dt 

d4K[Ste7ppFus3paut\ + p2if [Fus3ppo„tMAPif - Po^t] 

d[Fus3pa^tMAPK - Pout] ^ alK[Fus3poutmAPK - P^m] - (dlK + piK)[Fus3poutMAPK - P<,„t] 
dt 

d\Fus3poutSTe7pp] 

— ; = o4if Ste7pp][Fus3po„t - (d4jf + p4 jf ) Ste7ppFus3po„t 

dt 

d[FMs3ppo„t] ^ _a2K[Fus3ppout][MAPK - Pout] + d2K[Fus3pp„utMAPK - P„„t] + p4,K[Ste7ppFus3pout] + ** offK(C7 + C14 + C15 + 
dt 

G18 + C19 + G23 + C25 + G26 + C27 + B7 + B14 + B15 + S18 + S19 + S23 + B25 + B26 + B27) + fe49[FMs3ppi„] - fe50[FMs3ppo„t] 

d[Fus3ppoutMAPK - Pout] ^ a2K[Fus3ppout][MAPK - Pout] - (d2K + p2K)[Fus3ppoutMAPK - Pout] 
at 

d[MAPK - Pout] ^ _alK[Fus3pout][MAPK - Pout] + {dlK + plK)[Fus3poutMAPK - Pout] - a2K[Fus3ppout][MAPK - Pout] + (p^^k + 
dt 

d2.)[F„.3pp„„.MAP^ - --'1 + ^fSt^lI^]^ 
d\Fus3outSte7] 

— ^ = fe24[Ste7][F«s3o„t] - fe25[Fus3o„tSte7] 



= -Bl(onKKK[Stell] + * * ouk k k [Stellpp] + onKK[Ste7] + onjc [Fiis3o„t]) + offKKKBI +** offKKKBS + offKKBA +** 

at 

offKKBb + offKBG + ** oJ/kBT - fel6[G^i7Ste20]Bl + fel7Cl + ft22[5te5i„] - fe23Bl 

= onKKKBllStell] - offKKKB2 - onKK[Ste7]B2 - onK[Fus3out]B2 + offKKBS + ** offKKBU + offKB12 + ** offKB15 - 

at 

kie[Gi3jSte20]B2 + kl7C2 

= * * onKKK[Stellpp]Bl - onKK[Ste7]B3 - oriic [FMs3o„t]B3 + offKKB9 + ** oZ/kkBIO + offKBlS + ** offKBU - * * 

dt 

offKKKBS - /cl6[G;37Ste20]S3 + fcl7C3 
dB4 

— — = onKKBl[Ste7] - offKKBi - onKKK[Stell]Bi - * * onKKKB4,[Stellpp\ - onK[FusSout\Bi + of }kkkB& + * * of JkkkB^ + 
dt 

offKBie + ** offKBig - kl6[Gf!fSte20]BA + fel7C4 

— — = — * *onK K K[Stellpp]B5 — onKKKlStell]B5 — onK[Fv,s3out]B5 + * * of fKKKBW + of fKKKBll + offKB17 + * * offKBlS — * * 
dt 

of fKKBb - kl6[Gfi-tSte20\Bb + fcl7C5 

= - * *onKKK[Stellpp\B6 + onKBl{Fus3out\ - offKB& - onKKK[Stell\B& - onKK[Ste7\B& + offKKKB12 + * * offKKKB13 + 

dt 

offKKBW + ** offKKB17 - fcl6[G^7Ste20]S6 + fel7C6 

= _ * *onKKK[Stellpp\B7 - ouk K K[Sten]B7 - oriK KlSte7]B7 + * * of fK K K B14 + of fK K K B15 + * * of f'^ k BIS + of fK K B19 - * * 

dt 

offKB7 - fel6[G^75te20]B7 + fcl7C7 
dBS 

—— = -onK[Fus3a^t]B8-ioffKKK+offKK)B8 + onKKK[Stell]B4 + onKKlSte7]B2 + offKB20 + **offKB2-i-kl6[Gi3jSte20]B8 + kl7C8 
dt 
dB9 

—— = * * onKKKB4[Stenpp] - pkkB9 - (* * offKKK + offKK)B9 - onK[Fus3ov.t]B9 + onKK[Ste7]B3 + offKB21 + ** offKB25 - 
dt 

kl6lGp'rSte20]B9 + fel7G9 



dt 

* *onKKK[Stellpp]B5-{** offKKK + **offKK)B10-onK[Fus3out]B10 + offKB24 + **offKB27 + pKKB9-kl6[Gi}-fSte20]B10 + kl7C10 

^Ell = -{offKKK + ** offKK)Bll - onK[Fus3ov.t\B\l + onKKK[Stell]B5 + of fKB22 + ** offKB26 - fel6[Gfl7Ste20]Sll + fel7Cll 

dt 
dBl2 



dt 

-onKK[Ste7]Bl2 - (offKKK + offK)B12 + onKKK[Stell\Be, + ouk [Fus3o^t]B2 + offKKB20 + ** offKKB22 - fcl6[G;37Ste20]B12 + fel7C12 
dS13 



dt 

**onKKK[Stellpp]B&-{**offKKK + offK)Bl3-onKK[Ste7]Bl3+onK[Fu83a,,t]B3+offKKB2l + **offKKB2A-kl&[Gfi-iSte2Q]B13 + kl7Cl3 
= * * onKKK[Stellpp]B7 - (* * offKKK + ** offK)Bli - ouk K[Ste7\Bli + offKKB25 + ** off'j^^B27 - fel6[G^7Ste20]B14 + fel7C14 

^^^^ii = -{offKKK + ** offK)B15 - onKK[Ste7]Blb + on k k K[Stell]B7 + of f k k B23 + ** off'j^j^B2e - kie[G^-(Ste2a]Blb + fel7C15 
dt 

—— = - * *onKKKlStellpp]B16 - onKKKlStell]B16 - {offKK + of fK)BW + onKK[Ste7]BQ + onK[Fus3a^t]B4: + of fKKKB20 + * * 
dt 

offKKKB21 - fel6[G;37S'te20]B16 + fcl7C16 

= - * *onKKK[Stellpp]B17 - pkB17 - {* * offKK + offK)B17 - onKKK[Sten]B17 + onK[Fus3aut]B5 + offKKKB22 + ** 

dt 

offKKKB24 - fcl6[G^37S'ie20]B17 + fel7C17 
dB18 



dt 

-**onKKK[Stellpp]B18-{**off',^i^+**offK)B18-onKKK[Stell]B18+offKKKB26 + **offKKKB27+pKB17-kl6[G3')'Ste2a]B18 + kl7C18 
= - * *onKKKlStellpp]B19 - {offKK + ** offK)B19 - onKKK[Stell\B19 + onKK[Ste7]B7 + offKKKB23 + ** offKKKB25 - 

dt 

fcl6[G;37Ste20]B19 + kl7C19 
dB20 

= -{offKKK + offKK + offK)B20 + oriK K K[Stell]B16 + onKK[Ste7]B12 + onK[Fus3aut]B8 - /sl6[Ga7Ste20]B20 + fel7C20 

dt 

— ; — = »*onKKK[Stellpp]B16-{**offKKK+offKK+offK)B21-pKKB21+onKK[Ste7]B13+onK[Fus3o^t]B9-kl6[GB'rSte20]B21+kl7C21 

dt 
dB22 

—T- = -(offKKK + ** offKK +o//jc)B22 - {pk)B22 + onKKKlStell]B17 + onKlFus3a^,t]Bll - fel6[G^7Ste20]B22 + fcl7C22 
dt 

— ; — = -{offKKK + offKK + ** offK)B23 + ouk K K [Stell]B19 + onKK[Ste7]B15 - fel6[Grt7Ste20]B23 + fel7C23 
dt 

= **onKKKlStellpp]B17 -{**offKKK + **offKK + offK)B24-pKB24 + onK[Fus3aut]B10 + pKKB21-kl6[Gi3jSte20]B24 + kl7C24 

dt 

— ; — = * * onKKK[Stellpp]B19 - {* * offKKK + offKK + ** offK)B2b - pkkB25 + onKK[Ste7]B14 - fel6[G«7S'te20]B25 + kl7C25 
dt 

—r- = -{offKKK + * * of f'jijf + * * of fK)B26 + onKKK[Stell]B18 + pkB22 - fcl6[G^7Ste20]B26 + fel7C26 
dt 

— . — = * * onKKK[Stellpp\B18 - {* * offKKK + ** offKK + ** of fK)B27 + pkkB25 + pkB24 - klQ[G b1 Ste20]B27 + kl7C27 
dt 
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= -ClioriKKKlStell] + * * onKKKl3tellpp] + onKKlSte7] + onKlFus3out]) + offKKKC2 + **offKKKC3 + offKKC4 + ** 

at 

offKKCb + offKCe + ** offKCr + kl6[G0'rSte2O]Bl - fel7Cl 
dC2 

—— = onKKKCl[Stell] - of fKKKC2 - pkkkC2 - onKK[Ste7]C2 - onK[FusS„.^t]C2 + of JkkCS + * * of fKKCll + of fKC12 + * * 
at 

of fKC15 + fcl6[G^7Ste20]B2 - fel7C2 
d03 

—— = * « onKKK[Stellpp]Cl + pkkkC2 - onKK[Ste7\C3 - onK[Fus3ont]C3 + offKKCd + ** offKKCW + offKClS + ** offKCli - * * 
at 

offKKKd + fcl6[G;37Ste20]S3 - fcl7C3 

dCA 

—— = oriKKCllSter] - offKKd - onKKK[Stell]CA - * * onKKKCi[Stellpp] - onK[Fus3out]Ci + of fxKKCS + * * of fuKKCQ + 
at 

offKCie + ** offKCW + fcl6[G^7Ste20lB4 - fel7C4 
dG5 

— — = — * *onK K K[Ste\lpp]C5 — ouk K K[Ste\l]C5 — onK[Fus3out]Ch + * * of fKKKCW + of fKKKCll + offKClT + * * offKCli — * * 
at 

offKKCh + fel6[G,37Ste20]B5 - fcl7C5 

^ = - * *onKKK{Stellpp\C& + onjfCl[F«s3o„t] - offKCQ - onKKK[Stell\C& - onKK[Ste7]C& + offKKKC12 + * * offKKKC13 + 
at 

offKKCW + ** offKKCn + fcl6[G^7Ste20]B6 - fel7C6 
dC7 

—— = - * *onKKKlStellpp]C7 - on k K K[Ste\l]C7 - onKK[Ste7]C7 + ** offKKKCU + offKKKClb + ** offKKCW + offKKCW - * * 
dt 

offKC7 + kW[GpTSte2Q\B7 - kl7C7 
dCS _ 
dt ~ 

-pKKKCS-onK[Fus3out]C8-{offKKK+offKK)CS+onKKK[Sten]C4+onKK[Ste7]C2 + offKC20+**offKC23+kW[Gii-ySte20]BS-kl7C8 

^ = * * onKKKC4[Stellpp] - pkkC9 - (* * offKKK + offKK)C9 - onK[Fus3aut]C9 + onKK[Ste7]C3 + offKC21 + ** offKC25 + 
at 

PkkkCS + fcl6[G;37Ste20]B9 - fcl7C9 

= * * onKKK[Stenpp]C5 - {** offKKK + **offKK)CW-onK[Fus3out]CW + offKC24 + **offKC27 + pKKKCn+pKKC9 + 

dt 

fel6[G37Ste20]B10 - fcl7C10 

= -pKKKCll-(offKKK + **offKK)Cll-onKlFus3a^t]Cll + onKKKlStell]C5 + offKC22 + **offKC26 + kW[G0jSte2O]Bll-kl7Cll 

dC 12 

—T- = -PkkkCW - onKKlSLe7]C12 - (of fKKK + of fK)C12 + onKKK[Stell]CQ + onK[Fu83o^t]C2 + of fKKC20 + * * of fKKC22 + 
dt 

kW[Gi3iSte20]B12 - fel7C12 
dClS 

— — = * * ouKKK {Stellpp\C6 -{** offKKK + o//k)C13 - OUKK [Ste7jC13 + owit [Fus3„„f ]C3 + offKKC21 + ** offKKC24 + pkkkC12 + 

dt 

fel6[G37S'te20]B13 - fel7C13 
dG14 



dt 

**onKKKlStellpp]C7-(**offKKK + **offK)C14-onKK[Ste7]C14 + offKKC25 + **off'^^C27 + pKKKC15 + kl6[Gi3'rSte20]B14~kl7C14 
dCW 

= -pKKKC15-{offKKK + **offK)C15-onKK[Ste7]C15 + onKKKlStell]C7+offKKC23+**off',^,^C26 + kWlG^',Ste20]B15~kl7C15 

dt 
dCW 

= - » *onKKK[Stellpp]CW - onKKK[Stell]C 16 - (offKK + offK)Cl& + onKK[Ste7]C& + onK[Fu83a^t]C4 + of fKKKC20 + * * 
offKKKC21 + fel6[G37Ste20]B16 - fcl7C16 

= - * *onKKK{Stellpp\C17 -pkC17 -(**offKK + offK)C17 -onKKK[Stell\C17 + onK[Fus3out]Cb + offKKKC22 + ** 

dt 

offKKKC24 + kWlG0iSte2O]B17 - fcl7C17 
dCW 



dt 

-**onKKK[Stellpp]C18~{**off'jfj^+**offK)CW-onKKK[Stell]CW + offKKKC26 + **offKKKC27+pKC17+kl6[G^3'rSte2a]BW-kl7CW 
= - * *onKKKlStellpp]C19 - (offxK + * * o//jf )C19 - onKKK[Stell]C19 + onKK[Ste7]C7 + offKKKC23 + * * offKKKC25 + 

dt 

kW[GffiSte20]B19 - fel7C19 

^^^^ = -{offKKK + offKK + offK)C20-pKKKC20 + onKKK[Stell]CW + onKK[Ste7]C12 + onK[Fus3„ut\C8 + kW[G;3'lSte20]B20-kl7C20 

^^^1^ = * * onKKK[Stellpp\CW - (* * offKKK + offKK + offK)C21 - pkkC21 + onKK[Ste7]C13 + onK[Fus3out]C9 + pkkkC20 + 

fel6[G37Ste20]B21 - fel7C21 
dG22 

^— = -{offKKK + * * offKK + offK)C22 - (pKKK + Pk)C22 + OUK K K[Stell]C17 + onK[Fus3out\Cll + fel6[G^7Ste20]S22 - fel7C22 

dt 
dC23 

—T- = -(offKKK + o//kk + * * offK)C23 - PKKKC23 + onKKK[Stell]C19 + onKK[Ste7]C15 + kW[GpiSte20]B23 - fel7C23 

at 
dC24 



dt 

**onKKK[Stellpp]C17-(**offKKK+**offKK+offK)C24-pKC24+onK[Fus3aut]CW+pKKKC22+pKKC21 + kl6[Gf,jSte20]B24-kl7C24 
^2^5 = **onKKKlStellpp]CW-{**offKKK+offKK+**offK)C25-pKKC25 + onKK[Ste7]C14+pKKKC23 + kWlGffjSte20]B25-kl7C25 

= -{offKKK + * * o//^^ + * * offK)C26 - PKKKC26 + owK/fK [Stell]C18 + pjf C22 + kWlGpjSte20]B26 - fel7C26 

dt 
dC27 

—— = * * onif [StellppjClS - (* * offKKK + * * o//^^^ + * * offK)C27 + pic kkC26 + pkkC25 + pjc C24 + fcl6[G/375^te20]S27 - kl7C27 
at 



'^^^^^^'"^ = -fe22[Ste5i„] + fe23Sl 
dt 

'^^^dt^'"^ ^ -fe47[Fu53i„] + k48[Fus3out] + P^KlFusSpinMAPK - Pin] 

d[Fus3pi„] ^ p2f.[Fus3ppinMAPK - Pi„] - alKlFus3pi„]lMAPK - Pi„] + dlKlFus3pi„MAPK - Pi„] 
dt 

d[Fus3pinMAPK - Pir,] , „ , , , , 

-i , = alK[Fus3pi„]lMAPK - Pi„] - dlK[Fus3pinM APK - Pi„] - plKlFus3pi„MAPK - Pi„] 

dt 

d{Fus3ppi„] 

— jJ-i- — 1 = -k49[Fus3ppi„] + k5Q[Fus3ppout] - a2K[Fus3ppi„][MAPK - Pi„] + d2K[Fus3ppi„MAPK - Pi„] 

dj u.-,.,pi,,„MAJ A - P.„] ^ a2K\Fus3ppin][MAPK - Pi„] - d2K[Fus3ppiriMAPK - Pi„] - p2K[Fus3ppinMAPK - Pi„ 
dt 

d[MAPK - Pin] 



dt 

-alK[Fus3pin][MAPK - P,„] + {dlu + plK)[Fus3pinMAPK - Pin] - a2K[Fus3ppin][MAPK - Pin] + (p2k + d2K)[Fus3ppinMAPK - 

d[Stel2^^u..] k29.[Fus3pp,n] , 
Jt = k29 + [Fus3pp.n] - '=30[Stel2„e«„e] 

d\Fari\ fe33[5tel2„ctive]^ r , , , , 

= ^,,2 ro. To 12 - fc36 Farl F«s3ppi„ - fe35 Farl + fe37 Farlppi„ 

d\FcLT\'D'D' 1 

-i ^-i-^ = A;36[Forl][F«s3ppi„] - k37[F arlppin] - k38[Farlppin] + fe39[Farlpp„„t] - fe43[Farlppi„][Cdc28] + /c42[Farlppi„Cdc28] 

dt 

d[Cdc28] ^ fe42[i7'arlppi„Cdc28] - fe43[Farlppi„] [Cdc28] 

dt 

d\FarlppinCdc2f,] 

— — = -k42[Farlppi„Cdc28] + /c43 Far lppi„ Cdc28 

dt 

'^^^"''dr^""*^ " -fc4O[Farlppo„t][G07] + kil[Farlppa^tG /j-y] + fc38[Farlppi„] - k39[Farlppout] 

d[FarlppoutGi3i] ^ ^^Q[G0'y]lFarlpp„^t] - fc41[Farlpp„„tG^7] 
at 
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Abstract: The mating pathway in Saccharomyces cerevisiae is one of the best under- 
stood signal transduction pathways in eukaryotes. It transmits the mating signal from 
plasma membrane into the nucleus through the G-protein coupled receptor and the mitogen- 
activatcd protein kinase (MAPK) cascade. According to the current understandings of the 
mating pathway, we construct a system of ordinary differential equations to describe the pro- 
cess. Our model is consistent with a wide range of experiments, indicating that it captures 
some main characteristics of the signal transduction along the pathway. Investigation with 
the model reveals that the shuttling of the scaffold protein and the dephosphorylation of 
kinases involved in the MAPK cascade cooperate to regulate the response upon pheromone 
induction and to help preserving the fidelity of the mating signaling. We explored factors af- 
fecting the dose-response curves of this pathway and found that both negative feedback and 
concentrations of the proteins involved in the MAPK cascade play crucial role. Contrary to 
some other MAPK systems where signaling sensitivity is being amplified successively along 
the cascade, here the mating signal is transmitted through the cascade in an almost linear 
fashion. 
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INTRODUCTION 

Cells have to respond to changes in the environment and/or to the external stimuli. This 
is accomplished by signal transduction pathways which sense the signal, transduce it and 
induce necessary changes in the cell, e.g. in gene expression. One of the best understood 
signaling pathways in eukaryotes is the mating pathway in the budding yeast Saccharomyces 
cerevisiae (1,2). Extensive studies of the mating pheromone response have contributed much 
to the understanding of the mechanisms of several conservative biological modules (3), such 
as the G protein cycle (2,4) and the mitogen-activated protein kinase (MAPK) cascade (2,5). 
Genetic, biochemical and molecular analysis of the response have combined to establish 
basic principles of the signaling and regulation. Many important discoveries are made in the 
study of this pathway, for instance, the concept of a kinase-scaffold protein (6) and the role 
of regulator of G protein signaling (RGS) proteins in the pathway (7). 

The budding yeast can exist in either of the two types, MATa. or MATa. These two types 
of cells will mate when each one receives the mating signal, a peptide pheromone, secreted by 
the opposite type (a-factor by MATa. and a-factor by MATa). Once the pheromone binds 
to the seven-transmembrane-segment receptor in the plasma membrane (Ste2 in MATa. and 
Ste3 in MATa). the receptor is activated, which then activates the heterotrimeric G protein 
that couples to it (Fig. 1). The activated G protein transmits signal to multiple effectors, 
resulting in the beginning of the MAPK cascade, which is embedded in a scaffold protein 
Ste5. This cascade consists of three kinases: Stell (MAPKKK), Ste7 (MAPKK), and Fus3 
(MAPK). The activation of the cascade finally leads to the phosphorylation of Fus3. The 
phosphorylated MAPK then travels into the nucleus, and transmits the signal to downstream 
effectors, leading to preparation for mating, including the cell cycle arrest in Gl phase to 
assure synchronism of the mating partners, the induction of new gene expression necessary 
for mating, and the polarized growth in the direction of the pheromone source. 

Much qualitative and quantitative information in this pathway have been documented. 
With the increasing amount of experimental data and information, it is now possible to 
study this pathway quantitatively at a systems level. Several mathematical models have 
been employed to study this (8-12) and some other related systems (13-16), showing that 
mathematical modeling and simulation can be a powerful method in the analysis of functional 
and structural characteristics of biological pathways. 

We set up an Ordinary Differential Equation (ODE) model to describe the mating 
pathway in budding yeast. Although several models have been employed to illustrate 
mechanisms in the pathway, there has not been one that integrates all the known es- 
sential features with a comprehensive analysis of its dynamic properties. Some models 
were constrained to a single step (8,12), while others oversimplified the regulations and 
functions of the scaffold (10). In our model, biochemical interactions, induced gene 
expressions which feed back to the pathway, and translocations of key components such 
as the scaffold protein Ste5 are all considered. Results from our model are consistent 
with a wide range of experimental data. We then tested the current understanding of 
regulations of cellular responses and further explored the intrinsic mechanisms in the 
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pathway, with special interest in the role of the scaffold protein Ste5. We find that 
the shuttling of the scaffold and dephosphorylation of the MAP kinases cooperate to 
regulate the responses upon pheromone induction, and to help keeping the fidelity of 
the mating pathway. We further explored the mechanisms of the dose-response curves of 
this pathway, and elucidated the role of enzyme concentration. We found that instead 
of an ultrasensitive response as in some other MAPK cascade (17), the mating signal 
here is transmitted through the cascade in an almost linear manner due to negative feedback. 

THE MODEL 

We choose a mutant (TMYlOl), a MATa. type of Saccharomyces cerevisiae, as the main 
modeling subject. In this type of cells, the gene BARl is deleted. In a wild type MATa 
cell, the product of BARl can be excreted from the cells and cleave the a-factor. In order 
to simulate a continuous and constant a-factor treatment, we use this mutant in our model. 

The mating response can be divided into three modules in a temporal order: the activation 
of G protein cycle, the scaffold-depended MAPK cascade, and the downstream effects of 
activated MAPK (See Fig. 1). Viewing the response as a series of modules arranged in 
the temporal order can help to better understand the signaling process. Yet in our model, 
couplings and feedback between these modules are also taken into account. 

The activation of G protein cycle, a-f actor secreted by MATa binds to and hence 
activates the scven-transmembrane-segment receptor (Ste2) on the plasma membrane sur- 
face of MATa. Pheromone binding enhances monoubiquitination of the receptor, and the 
ubiquitination in this case serves as a signal for endocytosis and delivery to the vacuole (18). 
This comprises a negative feedback loop (at short time scales) . In our model, this process is 
treated as a process of accelerated degradation for simplicity. The synthesis of receptor Ste2 
is included; the downstream effecter Stel2 is responsible for the gene expression of Ste2. 
Thus, it comprises a positive feedback (at long time scales). 

The interaction between the activated receptor and Ga leads to some conformational 
changes, which enable Ga to release GDP and to bind GTP (19). Ga-GTP can not interact 
with G/37, resulting in a release of G(3'y from the receptor. The G7 unit fixes the heterodimer 
on the plasma membrane surface, while the G(3 unit can interact with several effectors to 
transmit the signal. In this sense, Ga unit is a negative regulator of the pathway; it plays 
a role in an adaptational response to pheromone through preventing the availability of Gf3^ 
when there is no signal (20). The Ga-GTP can be hydrolyzcd into Ga-GDP, which can 
re-associate with G(3j into a heterotrimer. The cycle of G protein is thus closed. Regulators 
of G protein signahng (RGS) proteins can accelerate the hydrolyzation of Ga-GTP (21). In 
this pathway, the most important RGS protein is Sst2, which is considered in our model. The 
gene expression of Sst2 is also regulated by Stel2 in the downstream. Therefore, Sst2 is part 
of a negative feedback loop that leads to the adaptation (8,22). Since there is experimental 
evidence that the amount of Ga increases significantly when the cells are treated with 
pheromone (8), we add G protein synthesis in our model. It is commonly accepted that 
Stel2 transcripts Ga genes. According to the above description, we formulate the reactions 
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in G protein cycle as follows: 



fel 



Ste2 + a — factor <^ Ste2acti- 



(degradation) 




(synthesis) 



G 




Ga ■ GTP + G/37 




fell 



Ga ■ GTP 



Ga ■ GDP 



Sst2<j(.tj„^,fcl4 



Ga ■ GDP + G/S-y ^ G 



where k's are the kinetic parameters, the protein above or below the arrow is enzyme or 
transcription factor of the reaction. 

The scaffold depended MAPK pathway. The released GPj has several effectors 
(23). One effector for mating is Ste20, the first p21-activated protein kinase to be identified 
in any eukaryote (24). Ste20 is also activated by Cdc42, which is regulated by Cdc24. 
However, this process is not included in our model because the Cdc42 binding domain of 
Ste20 has been shown to be dispensable for pheromone signaling in yeast (25,26), and there 
should be enough active Cdc24'^^^ and Cdc42 constitutively at the membrane to activate 
the amount of Ste20 required for initial signahng. Besides, mutants in Cdc24 do not have 
much influence on the pathway (27,28). 

Another effector of 0/3^ is the scaffold protein Ste5. The correlation between the disrup- 
tion of the Ste4(G/3)-Ste5 interaction and sterility confirms the importance of this interaction 
in signal transduction (29). G/?7 can bind to Ste5 on the LIM domain of Ste5, which is re- 
quired for Stell (MAPKKK) activation (30), probably through inducing a conformational 
change that enhances Ste20-dependent activation of Stell. Also it interacts with Ste5 in 
the R1NG-H2 domain which is essential for Ste5 oligemerization (31). 

Most scaffolds are contained in the nucleus during vegetative growth. Upon pheromone 
induction they undergo enhanced exportation from the nucleus and localize at the shmoo 
tip (1). Although the detailed controlling mechanism of exportation of Ste5 is not clear, it 
is plausible that mating pheromone increases the rate of Ste5 export (57). Here, we utilize 
an active control mechanism where the import rate is kept constant, while the export rate 
is dependent on the total concentration of the released Gf3'y. When there is no signal, the 
export rate is very low, keeping most scaffolds in the nucleus. When the mating signal opens 
the G-protein cycle, released G/S'-f enhances the export rate, driving scaffolds to the shmoo 
tip. In this way, the localization of the scaffolds can be regulated by G protein cycle. 

The mating pathway is highly dependent on the scaffold protein Ste5. First, Ste5 func- 
tions as an adapter protein. It recruits Stell to the plasma membrane, where Ste20 is 
also tethered, to facilitate Stell 's activation (28), triggering the MAPK cascade. Another 
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function for Ste5 is scaffolding. Stc5 tetliers Stell (MAPKKK), Ste7 (MAPKK) and Fus3 
(MAPK) to form a complex (32), keeping the kinases and their substrates in proximity, as 
well as preventing the influence of phosphatases. This function is supposed to be important 
in enzyme regulation and in preventing cross talk (33,34). 

When Ste5 is in the cytosol, it can form scaffold-kinase complexes with Stell, Ste7 and 
Fus3. Every kinase binding site on the scaffold is in one of the three possible states: without 
a kinase, with an unphosphorylated kinase, or with a dual phosphorylatcd kinase. So for 
scaffold-kinase complexes in solution, there are 3 x 3 x 3 = 27 states: Bl, B2, ...B27 (See 
Fig. 2A). G(3j can bind to Ste20 and Bi {i = 1,2. ..27). Because Ste20 is already on the 
plasmid membrane through the interaction with Cdc42 before signaling, and scaffold must 
shuttle out from the nucleus to bind to G/37, we assume that 0/^7 first binds to Ste20, then 
binds to Bi. Once Bi binds to Gf3'ySte20 complex, it is fixed at the plasmid membrane and 
the whole complex is denoted Ci (see Fig. 2B). Ci and Bi are the same in the interaction 
with MAPK kinases, except that Ci can phosphorylate Stell while Bi can not. 

fcl8 

G/S-y + Ste20 Gf3jSte20 
kl9 

klG 

Bi + Gl3-/Ste20 ^ Ci(i = 1, 2, ...27) 

fcl7 

fc22 

Ste5in ^ Bl 



where " Ste5in' denotes Ste5 in the nucleus, Bl denotes Ste5 outside of the nucleus, and k22 
is dependent on the total concentration of released G/37: A;22 = 0.0003 + 0.3 '^^2500 ' ^"ct — 
[G^-GTP] + [G^-GDP]. 

Stell is the MAPKKK of yeast pheromone pathway, which consists of an N-terminal 
regulatory domain and a C-terminal kinase region (1). The interaction of these two do- 
mains keeps Stell in an anti-self state. CBD domain in N-terminal contains serine and 
threonine residues that can be phosphorylatcd by Ste20. Ste20-mcdiatcd phosphorylation 
of these residues activates Stell (35). In addition to recruiting Stell to a pool of its ac- 
tivators Ste20, Ste5 also binds to Stell in its N-terminal, making the CBD domain in 
N-terminal more accessible for Ste20. SteSO also helps to make the CBD domain more ac- 
cessible to Ste20 by a direct interaction between the SAM domain of itself and the SAM 
domain of Stell (36), but is less essential than Ste5. Cells lacking SteSO are not truly sterile, 
thus we do not include SteSO in our model. Stellpp phosphorylates the target residues in 
the activation loops of Stc7 (MAPKK). and activates it (37). The activated Ste7pp then 
phosphorylates, and activates its targets, the MAPKs Fus3 and Kssl (38). In the mating 
pathway, Fus3 plays a much more important role, while Kssl is the main MAPK in the 
filamentation-invasion pathway in nitrogen starved cells (S). Thus we do not consider Kssl 
in our model, although analysis about crosstalk will be given in the Discussion. While Stell 
and Ste7 are predominantly cytoplasmic proteins, Fus3 can shuttle between the nucleus and 
the cytoplasm. It concentrates in the nucleus after activation, thus bringing the signal to 
the nucleus (7,3S,39). There are also several feedback loops: Stell (MAPKKK) undergoes 
ubiquitination and MAPK dependent degradation (40); Ste7 (MAPKK) is assumed to un- 
dergo enhanced degradation after phosphorylation (41). In addition, Ste7pp in the scaffold 
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is assumed to be hyper-phosphorylatcd by activated Fus3pp, which reduces the binding 
efficiency between Ste7pp and the scaffold sharply (42). 

All the kinases can be categorized into two pools: on the scaffold and in the solution. We 
assume that phosphorylation on scaffold employs a processive mechanism, while phospho- 
rylation in solution is distributive(ll). Processive mechanism means that the active kinase 
collides with and binds to a substrate, phosphorylates it once, then it may slides to align 
the second phosphorylation site of the substrate with the active site of the kinase, phos- 
phorylates the substrate a second time before finally dissociates. Distributive mechanism 
means that the active kinase collides with and binds to a substrate, phosphorylates it once 
and releases the monophosphrylated product, which then collides with a second molecule 
of the active kinase, and is phosphorylated a second time (14). In our model, we assume 
that a kinase is activated when and only when it is dual phosphorylated, while a partial 
phosphorylated kinase possesses no activity. Dual phosphorylation in a distributive manner 
could lead to a sharp, sigmodial stimulus-response curve (14,17), leading an all-or-none cell 
fate (14). However, the scaffold might diminish this property if phosphorylations on the 
scaffold occur in a processive manner (11). The dephosphorylations in the solution em- 
ploy the distributive mechanism while dephosphorylations in scaffolds are precluded in our 
model due to sterical obstruction of the phosphatases groups. The proteins responsible for 
dephosphorylation of Stellpp and Ste7pp are not clear. In the model, we add two proteins 
with constant concentration to dephosphorylate Stellpp and Ste7pp respectively. There are 
several phosphatases for Fus3pp: the dual-specificity phosphatase Msg5 (equally distributed 
in nucleus and cytoplasm), and the tyrosine phosphatases Ptp3 (cytoplasm) and Ptp2 (nu- 
cleus), all of which can result in the inactivation of FusSpp (7,43,44). The basal level of Fus3 
phosphorylation is controlled mainly by Ptp3, the amount of which is constant during the 
stimulation (44). Pheromonc treatment induces the expression of Msg5 through the effects 
of Stel2 (43), which then acts together with Ptp3 to inactivate Fus3pp. In our model, we 
use MAPK-P with an initial concentration and with a synthesis rate regulated by Stel2 
to represent these three phosphatases. A recent experiment shows that different inputs by 
Ste5 and Msg5 phosphatase lead MAPK cascade to multiple outcomes (45), indicating that 
MAPK-P is a key regulator in the network. Reactions of MAPK cascade in cytosol are 
formulated as follows: 

Stellp + MAPKKK - P ^ StellpMAPKKK - P -» Stell + MAPKKK - P 
Stellpp + MAPKKK -P ^ StellppM AP K K K - P -* Stellp + MAPKKK - P 
Ste7 + Stellpp ^ StelStellpp -* Stelp + Stellpp 

Ste7p + MAPKK - P ^ StelpMAPKK - P -* Ste7 + MAPKK - P 

Stelp + Stellpp ^ StelpStellpp -» Stelpp + Stellpp 

Stelpp + MAPKK ~ P ^ Ste7ppMAPKK - P ^ Ste7p + MAPKK - P 
Fus3ov,t + Ste7pp ^ Fus3outSte7pp — > FusSpout + Ste7pp 

FusSpout + MAPK - Pout ^ FusSpoutMAPK - Pout -* FusSout + MAPK - Pout 
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Fus3pout + Stelpp ^ Fus3poutSte7pp — > Fus3ppout + Stelpp 

FusSppout + MAPK - Pout ^ Fus^ppoutMAPK - Pout -» Fus^pout + MAPK - Pout 
/c24 

Stel + FusSout ^ Ste7Fus3out 

k25 



Ste7pp 

where p indicates once phosphorylation, and pp indicates twice phosphorylation. 

As for the scaffolds, wc made the following assumptions in the model: (a) Inactive kinase 
can bind to Bi and Ci. On the scaffold, this inactive kinase can either dissociate from the 
scaffold without phosphorylation or undergo processivc phosphorylation before getting off 
the scaffold if its upstream kinase happens to be on the same scaffold and in the active state, 
(b) Dephosphorylations on scaffolds are precluded due to sterical obstruction, (c) There is 
no binding of partially activated kinases to the scaffold proteins. For free fully activated 
kinases, only Stellpp can bind to the scaffold. Experimental evidence indicates that active 
Ste5 can also accept Stellpp activated by other pathways and channel those signals to Fus3 
(46), and that a greater amount of scaffold proteins interact with Stell rather than the 
other two kinases. As for Ste7, it undergoes hyper-phosphorylation by activated Fus3pp, 
which accelerates its dissociation from the scaffold (42). Fus3pp dissociates rapidly from the 
scaffold after phosphorylation (47) and travels into the nucleus. Thus, the reassociation of 
Ste7pp and FusSpp to scaffold seems to be unlikely, (d) Scaffold molecules do possess some 
catalytic properties (33), so that the reaction rates within a scaffold complex are greater 
than in the solution. Moreover, Ste7 and Fus3 can bind firmly (23), and the residues for 
Ste7 binding in Fus3 are the same as the residues for Ste5 binding (48) , so it is reasonable 
to assume that Ste7 competes with Ste5 for binding to Fus3. Fig. 3 illustrates the scaffold- 
dependent reactions of Stcf I. 

Downstream effects: After activation, Fus3pp dissociates rapidly from the scaffold, 
while the scaffold remains tethered to the plasma membrane and partly in the solution 
(47,49), acting as a platform for activation of many molecules of Fus3 and leading to the 
propagation of the signal. The activated Fus3pp transmits the signal into the nucleus, 
resulting in the activation of transcription and the induction of cell cycle arrest. FusSpp is 
assumed to mediate the pheromone-induced transcription of PRE-containing genes through 
phosphorylation and activation of at least three nuclear proteins: Digl, Dig2 and Stel2 
(1,2). In unstimulated cells, Digl and Dig2 bind to and thus repress Stel2 (50). Fus3pp 
phosphorylates Digl, Dig2 and Stcl2, and induces the release of Stel2 from the complex 
(50,51). The free Stcl2 then interacts with other proteins of the transcription machinery 
and thereby activates transcription of many different genes. Among the products of these 
genes are proteins that activate (e.g. Fus3, receptor) or inhibit (e.g. Msg5, Sst2) the 
pathway (2). Therefore, the transcription affords several feedback loops in the pathway. 
Another important substrate of FusSpp is Farl. Activated Stel2 increases the transcription 
of Farl, and FusSpp is able to phosphorylate Farl and thus to stabilize it (52). Farl is 
a bi-functional scaffold protein. In the cytoplasm, Farl is involved in polarized growth; 
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while in the nucleus, it has a key function in controlling cell cycle (53,54). Farl inhibits 
Cln-Cdc28 complex, the master regulator of the yeast cell cycle in Gl phase. In our model, 
the binding of Farl to Cln-Cdc28 is treated as a symbol for cell cycle arrest. Stel2 is also 
the transcription factor of Barl, which is then excreted from the cell and inactivates a- 
factor (53). It is deleted in our model because most experiment results that we use were 
from barlA mutants. However, barlA could cause more than 100-fold sensitivity increase 
in downstream transcription response (12). Reactions in the downstream are: 

Fus3ppi„,k28,k29 
k30 

Stel2^^U.e.k31,k6^ ^^^^ _ ^^^^ 
Stel2a,k32,k5 



FusSout 

Stel2<,^tj„5,fc33,fc34 Fus3ppi„,fc36 fc38 

^ Farl ^ Farlppin ^ Farlppout 

k35 k37 k39 

fc4() 

Farlppout + GfH-y ^ FarlppoutGji'y 

fc41 

fc42 

Farlppin + C'dc28 ^ Farlppi„Cdc2S 

fc43 

Stel2„(.ti„5,,fc44,fe45 
fc46 
fc47 

Fus3in ^ Fus3out 
48 

fe49 

Fus3ppi„ ^ FusZppout 

fcSO 

Fus3ppi„ + MAPK - Pin ^ FusSpp^nMAPK - Pi„ 
-» Fus3pi„ + MAPK - P,„ 

Fus3p^„ + MAPK - P,„ ^ Fus3p,nMAPK - Pin 
Fus3ir, + MAPK - Pi„ 



Spatial location: we consider two compartments in the cell: the nucleus and the shmoo 
tip (a projection towards the direction of pheromone formed as a result of polarized growth). 
The nucleus is where downstream effects take place. The shmoo tip is where the many 
signaling proteins are concentrated (55) and the main place for the upstream reactions, 
including the G-protein cycle and Ste5-related reactions. Thus in our model we neglected 
the rest part of the cytosol. In other words, we "restricted" the cytosol to the shmoo tip. 
The scaffold protein Ste5, the MAPK Fus3 (both activated and inactivated), and Farl are 
shuttling between the nucleus and the shmoo tip. Ste5 is mostly sequestered in the nucleus 
in the absence of pheromone while pheromone enhances nuclear exportation of Ste5 (56). 
Nuclear localization of Fus3 is slightly enhanced by the pheromone treatment (7,39,57,58). 

The mathematical model: we employ a set of ordinary differential equations (ODEs) 
to describe the changes in the concentration of proteins involved in the mating pathway. 
Generally, in a system of 1 biochemical species with the concentration q (i=l,2,- • ■ ,1) and m 
biochemical reactions with the rates vj (j=l,2,- ■ • ,m), the following series of equations can 
be utilized to describe the biochemical mechanism in the system: 

dci 

—r- = /l(ci,C2, • • • ,Ci) =nilVl+ ni2V2 H 1- "Im^'m 

dt 
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dC2 , , , 

— ;— = /2(C1,C2, • • ■ , C() = 7121^1 + 71221)2 H h n2mVni 

at 



dci . , , 

-TT = /i(ci,C2, •••,€()= nafl +ni2V2 H \-nimVrr 

at 



The quantity denotes the stoichiometric coefficient. The rate of a reaction is a func- 
tion of the concentrations of substrates, products and probable effectors (10). If we treat 
the gene expression as a special kind of reaction which can be described with Hill functions, 
the equations listed above can be employed to describe the dynamics of our system. In our 
model, all the unbound substances in various phosphorylation states and complexes formed 
by them are viewed as individual species. All complex formations, dissociations, degrada- 
tions, phosphorylations, and dephosphorylations are treated as reactions. The parameters 
and initial concentrations in the model are derived from experiments whenever possible. 
For the remaining parameters, some are determined by fitting the results of the model to 
indirect experiments; others are estimated according to the mechanisms and similar reac- 
tions in other organisms. The list of the model parameters as well as the detailed ODEs are 
presented in SUPPLEMENTAL MATERIAL. For simulation, we use Matlab, version 6.5 
(The MathWorks, Natick, MA). 

RESULTS 

A. Temporal characteristics 

G-protein cycle. Fig. 4 summarizes the dynamics of the G- protein cycle. Upon 
saturated pheromone induction (1 /iM a-f actor), the level of the activated G-protein climbs 
up rapidly, reaches its peak at about 30 s, and then gradually declines to a bottom at 
about 7.5 minutes before it gradually increases again, as shown in Fig. 4A. The simulation 
result (solid line) fits quite well with the experiment data (12) (circles with error-bars). 
One crucial factor that might contribute to enhancing the closure of G-protein cycle is the 
endocytosis of activated receptor Ste2. This hypotheses is supported by experiment with 
mutant Ste2'^^^^ (the C-terminal tail of the a-factor receptor gene STE2 is removed to impair 
its endocytosis) (12). We slowed down the degradation rate of the active Ste2 to simulate 
the Ste2^^^^ mutant. Gonsistent with the experiment data, the closure of G-protein cycle is 
apparently impaired and the amount of the activated G-protein levels off after reaching its 
peak, as shown in Fig. 4B. The behavior of Ste2^^^^ cells (dashed line for simulation, and 
up-triangle with error-bars for experimental data) indicates that endocytosis is a key factor 
that causes G protein cycle to close up. Fig. 4A shows that after about t=10 minutes, the 
activated G-protein continues to rise steadily. We attribute it to protein synthesis, because 
that is the time scale for gene expression. To test this hypothesis, we delete protein synthesis 
of all proteins considered in our model, and find that the level of the activated G-protein 
does not rise in the simulation. The behavior of cycloheximide treated cells, as shown in 
Fig. 4B (dotted line for simulation and square with error-bars for experiment data (12)), 
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supports this hypotheses. For comparison, time-course for TMYlOl cells is also shown in 
the figure (solid line for simulation, and circles with error-bars for experiment data). 

Binding of Ste20 to Gj3^: After activation, Qj3^ activates two effectors: Ste20 (MAP- 
KKKK) and Bi (Scaffold in the solution), hence transmitting the signal downwards. In our 
model, the time-course for the pheromone induced binding of Ste20 to Gj3'-f fits well with 
experimental data (24), (Fig. 4C). It shows that Ste20 binds quickly to G/57 during the first 
5 minutes. The binding slows down afterwards and then speeds up. This time-course seems 
to follow the activation of G- protein cycle upstream (Fig. 4A), consistent with the pre- 
sumption that G/37-dependent activation rather than Cdc42-dependent activation of Ste20 
is critical in the mating pathway. 

Activation of MAPK pathway: Following the recruitment of the scaffold protein 
Ste5 to the membrane, signal passes down through the MAPK cascade (Fig. 4D). Note that 
except for the gradual recruitment of Ste5, the signal transduction is very fast. 

Downstream effects: Fig. 4E shows the activation of Stel2, FarlppG/57 and 
FarlppCdc28 to illustrate the downstream effects of Fus3pp (MAPK). The formation of 
complex Farlpp-Cln-Cdc28 is responsible for the cell cycle arrest, and Farlpp-G/37 causes 
the polarized growth and the formation of the shmoo tip. Since G/37 is a part of the up- 
stream complex involving Ste5 to provide a scaffold for the MAPK pathway, excess G/37 is 
not available until the pathway is attenuated to some extent. Thus the curve for Farl-G/37 
begins to rise at 20 minutes after pheromone treatment, relatively late compared to other 
downstream effectors. 

B. Features of the pathway 

Scaffold shuttling and dephosphorylation cooperate to regulate MAPK cas- 
cade quantitatively and to keep its fidelity to mating signal: One of the most 
distinctive features of the mating pathway is its dependence on scaffold. Interestingly, we 
found that the amount of Ste5 localizing out of nuclear upon pheromone induction coincides 
with activated Fus3pp (MAPK) in the dose-response curve, as shown in Fig. 5A. Thus we 
speculate that the mating pathway is tightly controlled by scaffold protein Ste5 and its shut- 
thng. It is highly possible that with different concentrations of scaffolds out of nuclear, the 
efficiency of MAPK cascade varies. When scaffold concentration is relative low, the pathway 
efficiency should increase with scaffold concentration because of the catalytical function and 
spacial protection function of the scaffold. When scaffold concentration is too high, it may 
reduce the mating efficiency because effective complex suitable for signal transduction is 
hardly found (9). To investigate this possibility, we shut down the shuttling of the scaffold 
protein Ste5 and varied the concentration of total Ste5 at the shmoo tip from InM to 1000 
nM. For the levels of pheromone induction, we varied the concentration from 0.1 nM to 
1000 nM. Consistent with our expectation, the intensity of the output of MAPK cascade, 
indicated by the amount of the activated Fus3pp, first increases with scaffold concentration, 
then decrease, with the optimal scaffold concentration around 100 nM, as shown in Fig. 5B. 
Thus, within our simulation range, pathway efficiency monotonously increases with scaffold 



12 



available in the shmoo tip, and it is reasonable to suggest that scaffold proteins have the 
ability to quantitatively control the strength of signal transmission through the MAPK cas- 
cade. This function of Ste5 could be achieved by facilitating Stell's activation by Ste20 
through binding to Gf3'~f, concentrating MAPK cascade components, and sequestering them 
from inhibition by phosphatases (See below). 

Dephosphorylation is the other mechanism we speculate that might contribute to con- 
troUing the activation of the mating pathway. To test this hypothesis, we varied MAPK-P's 
concentration (both in the nucleus and at the shmoo tip) from 100 nM down to nM. As 
expected, Fus3pp (MAPK) shows super-sensitivity upon pheromone induction, see Fig. 5C. 
How does dephoshorylation control activation of MAPK cascade? Since the kinases are 
exposed to phosphatase only when they are in the cytosol and the scaffold could help to 
prevent the influence of the phosphatase on the kinases bound to it, we suggest that certain 
level of phosphatase concentration can keep the kinase phosphorylation in the cytosol at 
a very low level, and thus constrain the signal transduction on the scaffold. Hence, when 
the phosphatases are attenuated, a large amount of activated Stellpp, Ste7pp, and Fus3pp 
could be accumulated in the cytosol even at lower level of scaffold protein recruited to the 
shmoo tip, bypassing the control of scaffold protein Ste5. In short, the shutthng of the 
scaffold and the dephosphorylation of the MAP kinases cooperate to control the activation 
of MAPK cascade quantitatively. 

The cooperation between the scaffold shuttling and the MAPK dephosphorylation is also 
crucial to the speciflcity of the pathway. There are at least 5 MAPK signal transduction 
pathways in budding yeast (5), some of which share the same proteins, such as Stell and 
Ste7. A big puzzle is how specificity is achieved. To investigate the mating pathway's ability 
to isolate inappropriate signals leaking in from other pathways such as the filamentation- 
invasion pathway, we tested the behavior of some mutants. We set all Stell molecules in the 
dual phosphorylated state at t = and shut off MAPKKK-P to simulate the constitutive 
activation of Stellpp in invasive growth in the absence of mating pheromone. Prom the 
time-course curve of Fig. 6A, we observed that although Ste7pp (MAPKK) is activated to a 
relatively low extent, little Fus3pp (MAPK) is stimulated. We then set all Ste7 molecules in 
the dual phosphorylated state at t = and shut off MAPKK-P to simulate constitutive Ste7 
activation, and found that FusSpp (MAPK) could only be activated transiently in cytosol - 
the activation dropped down immediately (within ten seconds) (inset in Fig. 6B). This result 
is consistent with the experiments, revealing that persistent activation by constitutive Ste7pp 
fails to support FusS-dependent mating in the yeast (42). Further simulation by removing 
MAPK-P indicates that the later deactivation of Fus3pp is caused by dephosphorylation 
(Fig. 6B). Note that only the kinases in the solution are exposed to phosphatases, and 
the scaffold could help to shield from the infiuence of phosphatases. Thus, phosphatases 
constrain the signal on the scaffold. While on the scaffold, Ste7pp prefers to phosphorylate 
Fus3 instead of Kssl, the signal constrained on the scaffold will lead to the activation of 
Fus3, rather than Kssl, the main MAPK in the filamentation- invasion pathway in nitrogen- 
starved cells. That means that the activation of Fus3 strictly rehes on the scaffold protein. 
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When there is no pheromone induction, no scaffold protein is recruited to the shmoo tip, and 
the phosphatases inhibit the basal activation of Fus3. With pheromone treatment, the active 
scaffold proteins (the activation means recruitment to the plasma membrane in our model) 
help to assemble MAP Kinases cascade components and accelerate the activation of Fus3pp, 
the main MAPK in the mating pathway. When activation exceeds dephosphorylation, the 
signal will be transmitted into the nuclear. 

To conclude, when there is no pheromone induction, phosphatases repress the activation 
of the mating pathway and prevent inappropriate signals from leaking in. When pheromone 
exists, scaffolds are recruited to the shmoo tip by the activated G protein, gathering MAPK 
cascade components and sequestering them from phosphatases so that the mating signal 
can be transmitted downwards. Thus, mating pathway is highly depended on scaffolds. 
This conclusion is consistent with experiments (45). Although there are other factors that 
contribute to suppress the crosstalk between the two pathways (34), the mechanism outlined 
above could also play an important role. 

Desensitization to pheromone induction: The amount of activated Pus3pp 
(MAPK) decreases with time even when the cells are exposed to prolonged a- factor in- 
duction, as shown in Fig. 7A our wild type cell simulation (solid hue). This indicates 
desensitization effect. Desensitization is a key feature of the pathway, which enables cells to 
reenter the cell cycle to resume vegetative growth. We investigated the possible factors that 
might contribute to this desensitization, and found that multiple negative feedback loops, 
such as the degradation of Stell (MAPKKK) and Ste7pp (MAPKK), the synthesis of Msg5 
(MAPK phosphatase) and Sst2, should be the major cause. As shown in Fig. 7A, a wild 
type cell with all the negative feedback shows desensitization, but a mutant without these 
negative feedback does not (dashed line) : the activation of Fus3pp does not decrease even 
after one hour treatment with saturating pheromone. 

Another important cause is the negative regulation of Ste7's (MAPKK) binding abihty 
to scaffold by Fus3pp (MAPK). In the scaffold, Ste7pp, which undergoes feedback phospho- 
rylation by activated Fus3pp, dissociates more quickly from the scaffold (42), hence exposes 
itself to ubiquitination and degradation (41,59). This feedback can also accelerate the dis- 
assembling of scaffold complexes. In our model, we assume that the dissociation rate for 
Ste7pp on scaffolds with both Ste7pp and Fus3pp is larger than that for Ste7pp on scaffolds 
without Fus3pp. If we change former parameter to be the same as the latter one, activation 
of Fus3pp continues to rise up after prolonged stimulation, as shown in the dashed line in 
Fig. 7B. On the other hand, the dissociation rate for Ste7pp from the scaffold without 
Fus3pp has to be relatively slow to keep the intensity of the signal transduction. A mutant 
in which the dissociation rate for the normal phosphorylated Ste7pp (with no Fus3pp on 
the scaffold) is enhanced (dotted line in Fig. 7B) results in a very low intensity of the sig- 
nal transduction. Thus our model shows that differentiated binding abilities of Ste7 to the 
scaffold ensure the right behavior of the MAPK cascade. 

Sensitivity to pheromone: Aside from the temporal characteristics, the sensitivity 
to different levels of pheromone induction is another key feature of the signal transduction 
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pathway. We compared the dose-response curves predicted by our model to those observed 
in experiments (12,45), and there is a quantitative agreement (Fig. 8). While our model 
specifically simulates the ABarl strain, our results are also consistent with experiments of 
the wild type, taking into account the 100 fold sensitivity shift (8,60). In our study, cells are 
assumed to be treated with indicated concentration (O.OOlnM-1000 nM) of a-factor for 20 
min. The response of every component to a certain concentration of a-factor is represented 
by the maximum amount of that component. For the activated Ste2, the peak is obtained 
within seconds; for G protein activation (measured by the sum of Ga-GTP and Ga-GDP 
), the maximum amount is achieved at about 30 seconds; for Stellpp, Ste7pp and Fus3pp, 
the peak value appears at about 20 minutes. These different time scales are consistent 
with experimental observations and illustrate the characteristics of different components' 
activation. All these curves are normalized. We fit the dose-response curves of Ste2a, G 
protein activation, scaffold recruitment, Stellpp, Ste7pp and FusSpp with the Hill Function: 

A X [ml" 
[out] - ^ ^ 



where n is the Hill coefficient and th the threshold where the response reaches half its 
maximun (Table III). Note that the sensitivity to a-factor is well conserved throughout the 
whole pathway, from receptor Ste2 at the very beginning throughout the MAPK cascade. 

As stated before, our model separates the whole mating pathway into different modules. 
It is interesting to explore the dose-response curve of each module. The first module is ligand 

binding. The reason why the dose-response curve of the activated Ste2 is a Hill function 

fci 

with n ~ 1 is the reaction it takes. Consider the reaction 5*^62 + a ^ Ste2act, with the 

k2 

input concentration of a (on the left hand side) fixed at [a] and the total concentration of 
Ste2 fixed at [Ste2]o. At steady state: 

with Ai = [Ste2]o and thi = k2/kl, which takes the value thi = 5.0 with our choice of 
parameters kl and k2. The coefficients derived from our simulation with our whole model 
are n = 0.9, th = 9.1. The differences come from protein synthesis, degradation and signal 
dependent feedback. When these effects are deleted from the whole model, the simulation 
results agree perfectly with the analysis. 

We further explored the dose-response curve for G protein activation. Again, for sim- 
plicity, we do not take into account in the calculation the signal- dependent production, 
degradation and the feedback of the RGS protein (we treat the amount of Sst2 as a con- 
stant) and look for the steady state solution. The result of the calculation is: 



[G]o X [Ste2 



where kl3' = A;13-|- kl4:[Sst2] . To get this result, we made the assumption that the 

fcl3' 
fcl5 



hydrolysis is a relatively slow process, so that << [G]o- If we substitute Eq.(l) into Eq. 
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(2), we get: 

where A2 = [G]oT^feiF ^^2 = thi tx^^tift- Eq. (3) is a Hill function with n = 1.0 

and t/i = 2.9. The simulation result of our whole model is n = 1.0 and th = 4.8. Again, 
if the protein synthesis, degradation and feedbacks are deleted, the simulation results of 
n = 1.0 and th = 2.9 agree well with the analysis. 

We see that the curve of [GP'y] ^ [a — factor] employs the same type of function as the 
curve of [Ste2act] ~ [a-factor], with the same Hill coefficient. The only difference is that 
there is a shift in threshold. 

The above analysis suggests that the amount of Sst2 can affect the threshold of the dose- 
response curve. Fig. 9A indicates that Sst2, which can accelerate the closure of the G 
protein cycle, is indeed a key regulator of mating pathway's sensitivity. Comparison of the 
dose-response curves for mutants SST2A (dotted line for simulation, up-triangle with error- 
bars for experiment data), wild type cells (solid line for simulation, squares for experiment 
data (8)) and 2x SST2 (dashed line for simulation, down-triangle for experiment data) 
clearly shows that the system is sensitive to the amount of Sst2, which is consistent with the 
theoretical analysis above. Another mutant we studied in this module is one with excess G(3j 
copies. Fig. 9B compares dose-response curves for TMYlOl cells (sohd hue for simulation, 
and squares for experiment data) and cells with 2x G/37 (dashed line for simulation, and 
up-triangles for experiment data(8)), indicating that G(3^ alone is sufficient to switch on the 
downstream signal transduction. 

The next module in the pathway is the scaffold- dependent phosphorylation cascade. 
Again, we do not take into account protein synthesis and degradation, nor all the feed- 
back. Here, GPj is the input, while the concentration of Stellpp, Ste7pp and Fus3pp is 
chosen as the output. 0/^7 first binds to Ste20, then binds to scaffold in the solution to 
generate C (we use C to indicate the ensemble of CI, C2,...C27) which then initiates MAPK 
cascade. We calculate the Hill coefficients of Stellpp, Ste7pp, and FusSpp using G/37 as 
input, and find that each of them shows ultrasensitivity (Numerical results: risteiipp — 2.3, 
i^steVpp = 2.2, njus^pp = 1.8). Since Gf3'~f experiences two simple reversible reactions to gen- 
erate C, the relation between C and Gf3^ must be a hill function with n = 1. Therefore, 
the ultrasensitivity in this module must arise from C to Stellpp, Ste7pp and Fus3pp. Let 
us first consider the phosphorylation and dephosphorylation cycle of Stell. Stell can be 
dual phosphorylated on the scaffold in a processive mechanism and dual dephosphorylated 
in the cytosol in a distributive mechanism: 

C + Stell ^ C Stell ^ C + Stellpp. 

off 

stellpp + MAPKKK - P ^ StellppMAPKKK -P^ Stellp + MAPKKK - P 
Stellp + MAPKKK - P ^ StellpMAPKKK - P ^ Stell + MAPKKK - P 
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where we assume that Stellpp leaves the scaffold so fast that CStellpp can be neglected in 
the process. Introduce Michaelis constant K = 2n±p^ Ki = K2 = and define 



ai ' ^ a2 

ai = eK± a2 = p' = K' = EiI<2±P2l£i we get: 



[Stellpp] = -^G (pCtotaup'[MAPKKK - P]o, tw^. jJ^) [Sten]o, (4) 



where 

2uN 



G {u, V, M, A^) 



V - u + vM + uN + ^{v -u + vM + uNf -A{v -u)uN 



u = p[^tell]o, V = p'[MAPKKK - P]o, M = ^ , ^ 



[5tell]o [5tell]o 
The above G function can be fitted to a sigmoidal curve (Hill function) (61), with the Hill 
coefficient and the threshold value to be 

1 



logs 



81(M+0.1)(Af+0.1) 
=81 (M+0.9)(iV+0.9) 



According to the parameters in the Stell cycle, M = K/[Stell]Q ~ 0.125 and A^ = 
K' /[Stell]o = 0.255. So Eq. (7) approximates a Hill function with n f» 2.6 and th ^ 21nM. 
Therefore, it is the zeroth-order ultrasensitivity that leads to the ultrasensitive response 
of Stellpp. Similarly, in the Ste7 cycle, M ^ 0.183 and A^ = 0.324; in the Fus3 cycle, 
M ^ 0.135 and A^ = 0.181. Therefore, Ste7 and Fus3 are also located in the zeroth-order 
region, making them ultrasensitive. 

The analysis above together with Eq. (6) indicates that the relation between Stellpp, 
Ste7pp, Fus3pp and a should exhibit ultrasensitivity (numerical result: Usteiipp = 2.1, 
nsteipp = 1-8, UfusSpp = 1-9). This appears to contradict the results in Table HI where all 
the dose-response curves overlap with Hill coefficient n ^ 1. The reason for this contra- 
diction is that in the analysis we cut off all the feedback in the original model. There 
are all together nine feedback in the whole model, six of which are negative feedback: the 
transcription of Sst2, MAPK-P, the degradation of activated Ste2, Stell, Ste7pp, and the 
hyper-phosphorylation of Ste7pp. When [a] is low, all these negative feedback are kept 
low, leaving the output nearly unaffected. When [a] is high, the negative feedback will 
also be strong, significantly reducing the output. Therefore, negative feedback can make the 
dose-response curve less steep. It may be the counterbalance between the zeroth-order ultra- 
sensitivity and the negative feedback that keeps the Hill coefficients of Stellpp, Ste7pp, and 
Fus3pp remain 1. In order to test this view, we add each of the six negative feedback into 
the simplified model where all feedback are cut off (Table |H] Numerical result shows that 
each feedback can reduce the Hill coefficients, and if all the negative feedback are added 
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together, the Hill coefficients can decrease to nearly 1, indicating that negative feedback 
indeed can reduce the Hill coefficients. 

We further study how the Hill coefficients of Stellpp, Ste7pp, Fus3pp depend on the scaf- 
fold, the substrate and the phosphatase concentrations (Table UTTl First, we increase Stell, 
Ste7 or Fus3 concentration by 10 times, and find that the Hill coefficients only change a little. 
The probable reason is that the original kinase concentration is already much larger than 
the scaffold concentration, so when Stell, Ste7 or Fus3 concentration increases, the added 
part could not get to the scaffold to be phosphorylated, and thus does not contribute to 
response. Then, we increase kinase concentration and scaffold concentration together, and 
find that Hill coefficients have a substantial increase. This is because when scaffold and ki- 
nase concentration increase together, effective substrate concentration increases. Therefore, 
the zeroth-order ultrasensitivity becomes more significant, while feedback, which is mainly 
dependent on the downstream regulation, does not increase as fast. Table UTTl also shows that 
when each phosphatase concentration is decreased by 10 times. Hill coefficients arise as well. 
This is because low phosphatase concentration allow more substrate to be phosphorylated 
in the cytosol in a distributive mechanism, which is regarded as another mechanism to gen- 
erate ultrasensitivity aside from the zeroth-order ultrasensitivity (17). Thus, the scaffold, 
the substrate and the phosphatase concentration play an important role in determining the 
Hill coefficients of the MAPK pathway. 

Analysis of parameters: Due to the lack of experimental data to determine all the 
parameters, it is necessary to analyze the sensitivity of the system to changes of the param- 
eters. To do this, we define a quantity: 



where [Fus3pp]cai denotes the calculated output when a parameter is changed and 
[Fus3pp]obs denotes the original output, and the bar denotes the average of the relative 
variance of the output over an input (a-factor) range of 10-3 nM to 10^ nM. We multiply 
and divide one parameter by 2 at each time, calculate D, and then take the average for the 
two situations of increasing and decreasing the parameter. The most sensitive parameters 
are listed in Table IIVI Note that these parameters all have direct or indirect experimental 
support, which is reassuring. These parameters also give some clues about the mechanisms 
of the mating pathway. The parameters which influence the outcome of the pathway most 
are those involved in the receptor activation, which is consistent with our finding that the 
shape of dose-response curve is determined by the ffist step in the pathway-the receptor ac- 
tivation. The parameter in the production of receptor Ste2 is also essential to the outcome 
of the system. The other two influential parameters are the rates of the kinases to get off 
the scaffold, so the scaffold protein is also an important factor in this pathway. 

DISCUSSION 

Our model describes the whole mating pathway comprising of the G protein cycle, the 
scaffold-dependent MAPK cascade and the downstream effects in the nucleus. We have in- 
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vestigated multiple features of the pathway, including its various characteristic times scales, 
dcscnsitization, scaffold's effect, specificity, sensitivity to different levels of pheromone in- 
duction, the role of feedback, and sensitivity amplification. Although many the parameters 
in our model do not have solid experimental support and the detailed mechanisms of some 
steps are still not clear, the results given by our model are consistent with the current 
understanding of the pathway and with a wide range of experimental data. 

The duration and sensitivity of the mating pathway have to be tightly regulated; an 
inappropriate activation of Pus3 will block the normal invasive growth. Our model shows 
that activation of the mating pathway attenuates with time even when the pathway is 
exposed to prolonged pheromone induction. This dcscnsitization is attributed to several 
feedback such as the enhanced degradation of Ste7 (MAPKK). These feedback enable the 
cell to recover from mating and continue their vegetative growth upon prolonged pheromone 
induction. 

Evidence shows that oligomerization of the scaffold protein is required for its activation 
(31,57,62,63). However, oligomerization is not included in our model for simplicity. Since we 
consider all possible complexes involving the scaffold protein, taking the oligomorization into 
account would make the model extremely complicated (with 27x27=729 possible complexes 
involving the scaffold protein). Moreover, there are few experimental data concerning the 
interactions between the scaffold protein and the kinases. Therefore, we have to make some 
simplifications. A clue for the simplification is the experimental evidence that nuclear export 
and shmoo tip recruitment of Ste5 are coordinated with oligomarization (57,61). Thus we 
use nuclear export as a controlling step. In our model, the recruitment of the scaffold protein 
to the shmoo tip implies its activation, including the effect of nuclear transportation and 
oligomorization. However, oligomarization may have other effect in the pathway besides the 
activation of scaffold protein. Due to the lack of experimental data in this process, it is 
difficult to consider it in detail in our model. Experimental progress about this process is 
much needed to further improve the model. 

The scaffold protein undergoes continuous shuttling and enhanced exportation upon 
pheromone induction. The obvious question is why yeast cells take so much trouble to 
shuttle a huge protein through the nucleus when it functions predominantly in cytoplasm. 
Our model indicates that nuclear shuttling might be a key step controlling the availability of 
the scaffold protein to the pathway. Since the activation of the MAPK cascade in the mating 
pathway is dependent on the scaffold protein, whether and how many scaffold proteins are 
available determine whether and how efficient the MAPK pathway is stimulated. However, 
the mechanism of scaffold shuttling are still not clear. The complex Msn5p/Ste21 is sug- 
gested to be responsible for the export of Ste5p. Further work is required to establish the 
accurate and detailed mechanism of this controlling step. In our model, we employ an active 
control mechanism. That is, the scaffold export rate is dependent on the concentration of 
the separated Gf3'~f, which is released by the phemorone. The more the pheromone, the more 
the activated GPj, and the higher the export rate. In other words, G protein cycle controls 
the shuttling of the scaffold and the concentration of scaffold at the shmoo tip. 
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Another function of the scaffold is its role in keeping the pathway's fidelity to the signal. 
The capability of the scaffold protein to prevent kinases from dephosphorylation assures the 
mating pathway's dependence on the scaffold protein, and the availability of scaffold in the 
shmoo tip is further controlled by the G protein cycle. Recently the specificity of different 
pathways in yeast is under intensive study (34). It is interesting to see that different cellular 
signals, which can be transmitted by the same components, result in distinct responses. Es- 
pecially, the haploid invasive growth pathway employs the same MAPK components (Stell 
as MAPKKK, Ste7 as MAPKK, Pus3 and Kssl as MAPK) as the mating pathway, except 
that FusS is more active during mating while Kssl is preferentially activated during invasive 
growth. Then how are the different outputs controlled? Our model suggests that dephos- 
phorylation and scaffolds work in coordination to prevent improper signal from leaking in, 
and thus contributing to the mating pathway's fidelity to pheromone induction. Further 
work is needed to include the parallel pathway of the invasive growth, through which a more 
comprehensive understanding of specificity might be obtained. 

MAPK cascade which is conserved in all eukaryotic cells is composed of both phosphory- 
lation and dephosphorylation. Our model reveals that dephosphorylation has several roles 
in the mating pathway. It is obvious that it contributes to the desensitization of the pathway 
which enables the cell to reenter the cell cycle. Furthermore, it cooperates with the shuttling 
of the scaffold proteins, to realize other important features of the signaling pathway. First, 
it helps to preserve the consistence of the sensitivity from the G-protein cycle to the MAPK 
cascade. Second, the amount of MAPK-P, together with the scaffold protein, contributes to 
the pathway fidelity. 

Notably, we suggest that negative feedback plays an important role in the experimentally 
observed preservation of sensitivity along the MAPK cascade. Whether or not the sensitiv- 
ity is amplified as the MAPK cascade descends is also determined by the concentrations of 
the kinases and phosphatases involved. We call for new experiments to test this hypothesis. 
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TABLE I: Coefficients in Hill function 



Coefficient Ste2a G activation Stellpp SteTpp FusSpp 

n 0.9 1.0 1.2 1.2 1.2 

th (nM) 9.1 4.8 4.0 4.8 4.8 



Note: simulation results with the whole model. 



TABLE IL Hill coefficient n when negative feedback are added 



added feedback 


Stellpp 


SteTpp 


FusSpp 


none 


2.1 


1.8 


1.9 


transcription of Sst2 


1.8 


1.5 


l.T 


transcription of MAPK-P 


2.0 


l.T 


1.8 


degradation of Ste2act 


1.8 


1.4 


1.8 


degradation of Stell 


1.3 


1.2 


1.8 


degradation of Ste7pp 


2.1 


l.T 


l.T 


hyper-phosphorylation of SteTpp 


1.8 


1.6 


1.4 


all the six feedback 


1.2 


1.2 


1.2 



TABLE III: Hill Coefficient n when concentration changes 



concentration changes 


Stellpp 


Ste7pp 


Fus3pp 


[Stell] X 10 


1.2 


1.2 


1.3 


[SteJ] X 10 


1.3 


1.6 


1.2 


[Fus3] X 10 


1.2 


1.2 


1.2 


[Stell] X 10, [Ste5] x 10 


1.7 


1.8 


1.4 


[Ste7] x 10, [Ste5] x 10 


1.8 


2.5 


2.8 


[Fus3] X 10, [Ste5] x 10 


1.4 


1.7 


1.7 


[MAPKKK - P]/W 


1.7 


1.7 


1.4 


[MAPKK - P]/10 


1.2 


1.7 


1.2 


[MAPK - P]/10 


1.2 


1.3 


2.1 



TABLE IV: Influential Parameters 



Parameter 

kl 

k2 

k6 

* * of Ik 

* * offxR 



D = 

0.519 

0.266 

0.167 

0.156 

0.150 



/ A[Fms3pp] \ 2 



Related reaction 

kl 

Ste2+a ^ Ste2act 

k2 



k4:,k5,stel2a ^, - fe7 
^ btez > 



fc6 



Fus3pp gets off from the scaffold protein 
Ste7pp gets off from the scaffold protein 
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FIGURE LEGENDS 

Figure 1. Spacial structure of the mating pathway. 

Figure 2. A. 27 solution-located scaffold complexes. B. 27 membrane-located scaffold 
complexes. The diamond on the upleft corner indicates the G/37Ste20 complex 

Figure 3. Scaffold-dependent reactions. 

Figure 4. Time course of the G protein cycle activation. A. G protein activation. The 
values are normalized to the maximum concentration around 30 s. Result from simulation 
is shown in solid line and experiment data (12) arc plotted in circles with error-bars. 
B. G protein activation in Ste2^^^^ cells, the wild types and cycloheximide treated cells 
(experiment data are from (12)) . C. The time-course of Binding Ste20 to G/37. The values 
are normalized to the maximum concentration (experiment data are from (24)). D. The 
recruitment of scaffold protein Ste5 (circles) which is the sum of scaffolds in the solution 
and scaffolds at the membrane, the activation of MAPK cascade components Stellpp 
(MAPKKKK) (dotted line), Ste7pp (MAPKK) (dashed line) and Fus3pp (MAPK)(sohd 
line). E. Downstream responses to a — factor induction: activated Stel2 (solid line), 
Farlpp ■ Cdc28 (dashed line), and Farlpp ■ G(3j. 

Figure 5. A. Predicted dose-response curves of recruitment of scaffold protein Ste5 
(dashed line) and activation of Fus3 (MAPK) (sohd hue). B. The dependence of Fus3pp 
on the concentration of scaffold proteins Ste5, with different concentration of a-f actor. C. 
Predicted dose-response curves for mutants in which MAPK-P is under-expressed. All 
the values arc normalized. In all above simulation. Cells are treated with indicated with 
a — factor for 20 min. 

Figure 6. Dcphosphorylation prevents improper signal to leak in. A. All Stell is dual 
phosphorylated at t = 0, and MAPKKK-P is shut off. This simulates the condition in 
which signal in invasive growth pathway is on. B. All Stell is dual phosphorylated at 
t — 0, and MAPKK-P is shut off. The activation of Fus3 is still repressed down except for 
the small pulse at the very beginning (Inset graph). Dashed hue indicates the activation of 
Fus3 when MAPK-P is eliminated. 

Figure 7. Descnsitization. A. Time-course of activation of Fus3 (MAPK) in a wild type 
cell (solid line) and in a mutant (dashed line). B. Effect of feedback hyper-phosphorylation 
of Ste7 in the scaffold. Predicted time-course activation of Fus3 for wild type (solid line) 
and for mutants (dashed line and dotted line). 

Figure 8. Dose-response curves for key components in mating pathway in TMYlOl cells: 
activated Ste2 (dotted line for simulation, down-triangles for experiment data from (12)); 
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G protein activation, i.e. the sum of Ga-GTP and Ga-GDP (dashed hne for simulation, 
circles for experiment data from (12)); Fus3pp (MAPK) (solid line for simulation and 
crosses for experiment data from (45)). 

Figure 9. Response curve of Stel2 to a-factor for: A. different Sst2 expression level; 
B. different G/37 expression level. To simulate the 2 x Gp^ cells, we separately add 
another lOOOnM G/37 at t=0 min to double the total concentration of G/37. Cells 
are treated with a-factor for 60 min. Note that the experiment data are shifted left by 
more than 100- fold because Barl is deleted in TMYlOl cells. Experiment data are from (8). 
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